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--H . ABSTRACT 

^ ■ 

' I As part of our radial velocity planet-search survey performed with SARG at TNG, we monitored the components of 

l/^ . HD 132563 for ten years. It is a binary system formed by two rather similar solar type stars with a projected separation 

of 4.1 arcsec, which corresponds to 400 AU at the distance of 96 pc. The two components are moderately metal-poor 

I I' ([Fe/H]=-0.19), and the age of the system is about 5 Gyr. We detected RV variations of HD 132563B with period 

Ph I of 1544 days and semi-amplitude of 26 m/s. From the star characteristics and line profile measurements, we infer 

Py , their Keplerian origin. Therefore HD 132563B turns out to host a planet with a projected mass msini = 1.49 Mj 

at 2.6 AU with a moderately eccentric orbit (e = 0.22). The planet around HD 132563B is one of the few that are 

known in triple stellar systems, as we found that the primary HD 132563A is itself a spectroscopic binary with a 

period longer than 15 years and an eccentricity higher than 0.65. The spectroscopic component was not detected in 

Q ■ adaptive-optics images taken with AdOpt@TNG, since it expected at a projected separation that was smaller than 0.2 

$_( ■ arcsec at the time of our observations. A small excess in K band difference between the components with respect to 

"ti ' the difference in V band is compatible with a companion of about 0.55 Mq. A preliminary statistical analysis of when 

j^ I planets occur in triple systems indicate a similar frequency of planets around the isolated component in a triple system, 

I I , components of wide binaries and single stars. There is no significant iron abundance difference between the components 

(A[Fe/H]=0.012 ± 0.013 dex). The lack of stars in binary systems and open clusters showing strong enhancements 

^"^ ' of iron abundance, which are comparable to the typical metallicity difference between stars with and without giant 

^ I planets, agrees with the idea that accretion of planetary material producing iron abundance anomalies over 0.1 dex is 

OO , rare. 

T-H . 

Qs , Key words. (Stars:) individual: HD 132563 - Planetary systems - (Stars:) binaries: visual - (Stars:) binaries: spectroscopic 

(— ^ ■ - Techniques: spectroscopic 



1. Introduction into the impact of this on the planet formation mechanism 

7—{ , , can be derived (Malmberg et al. I2007P . 

• • , Planets in multiple stellar systems play a special role in . r i ■ • ■ ..... 

> ■ determinations of the frequency of planets in the Galaxy, , ^^ ^P^^^ °f selection criteria against binaries m many 

k> ■ because more than half of solar type stars are in multiple planet-search surveys, several tens of planets are current 

^ : systems (Duquennoy & Mayor [IMID. Furthermore, plan- ^nown m mu tiple syst«ns (Desidera & Barbieri WM 

C^ ■ ets in binaries allow one to better understand the processes Mugrauer & Neuhauser mm. A variety of binary config- 

of planet formation and evolution, thanks to the dynami- orations are observed: several of the planets are m rather 

cal effects induced by the stellar companions on the proto- ^^^^^ ^'""f ""'' °' ^ "" *^ ''''''*' caused by the compamon 

planetary disks and on an already formed planetary system ^^ ^^^^ planetary system are expected to be small. However, 

(Marzari et al. WM • Finally, since the frequency of (wide) t^ereare few planet hosts with companionsal a separation 

binaries seems related to the density of the formation envi- t^^tissmall enough (e^g. Hatzes et al WM Chauvm et 

ronment (Lada & Ladal^Dmi Goodwin [MS some insight ^^^ ^^ *« ^^^^^ i* difficult to allow planet formationac- 

cordmg to our current knowledge (see e.g. Thebault I201ip 

even if they are on dynamically stable orbits. The variety of 

Send offprint requests to: S. Desidera, system configurations is further extended by the discovery 

e-mail: silvano . desideraSoapd . inaf . it of planets in triple systems and, more recently, of planets in 

* Based on observations made with the Italian Telescopio circumbinary orbits around close eclipsing binary systems 

Nazionale Galileo TNG) operated on the island of La Palma by ( j^ee et al-Wm, when using the timing technique. 

the Fundacion Galileo Galilei of the INAF (Istituto Nazionale ^ ^ ^' ^ 6 h 

di Astrofisica) at the Spanish Observatorio del Roque de los Planet properties also appear to be affected by the pres- 

Muchachos of the Instituto de Astrofisica de Canarias. ence of companions, as is the case for the mass distribution. 
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with massive planets in close orbits found primarily in bi- 
naries (Zuckcr & Mazeh l 2002l Udry et al. l2003l Desidera 
& Barbieri 2007), and for planet eccentricity, which is en- 
hanced in binary systems, possibly through the Kozai reso- 
nance (Tamuz et al. I2008P . The circumbinary planets found 
so far around eclipsing binaries are all rather massive and 
have long periods. The possible link of these features with 
the original system configuration and its evolution (all the 
central binaries include a remnant object) needs to be in- 
vestigated. 

The general radial velocity (RV) surveys include several 
binaries in their samples, but in most cases only the bright- 
est star in the system is under RV monitoring. Dedicated 
surveys of planets in binary systems with well-defined prop- 
erties help for better understanding the impact of binarity 
on planet formation. Some of these projects are focused 
on spectroscopic binaries (Eggenberger 120101 Konacki et 
al. 1^09)) , while others study visual binaries using RV (our 
survey, Desidera et al. 120101 Toyota et al. l2009p and astrom- 
etry (Muterspaugh et al. I2010p . Our ongoing survey using 
SARG at TNG is the first planet-search survey specifically 
devoted to binary systems. It targets moderately wide bi- 
naries (typical separation 100-500 AU) with similar compo- 
nents (twins) . The sample includes about 50 pairs of main 
sequence solar-type stars with projected separation larger 
than 2 arcsec, distance smaller than 100 pc and magnitude 
difference AV < 1 mag. The sample selection and the sur- 
vey itself are described in Desidera et al. (|2010p . 

Unlike most of the multiple systems included in the 
RV surveys, both components of the pair are systemati- 
cally observed in our survey, usually with the same time 
sampling. This allows the investigation of possible factors 
driving planets in binary systems to occur and of possible 
effects induced, directly or indirectly, by the presence of 
planets, such as altering chemical abundances by accreting 
metal- rich planetary material (Desidera et al. I2004|) . 

In this paper we present the interesting case of the HD 
132563 system, whose components show both RV varia- 
tions: those with amplitudes larger than 1 km/s for the pri- 
mary, implying most likely a stellar companion, and those 
with amplitudes of a few tens of m/s for the secondary, 
suggesting there is a planetary companion. The paper is 
organized as follows. In Sect. [5] we present the observations 
taken at TNG using SARG (high-resolution spectroscopy) 
and AdOpt (direct imaging). In Sect. |3] we present the prop- 
erties of the components. In Sect.|3]we present the evidence 
for a planetary companion around HD132563B. In Sect. [5] 
we show that HD132563A is a long-period spectroscopic 
binary, and we discuss the clues to the mass and the or- 
bital parameters of the companion. In Sect. [5] we present 
the clues to the orbital parameters of the wide pair and 
its possible infiuence on the planet orbiting HD 132563B. 
In Sect. [7] we discuss the results, with emphasis on planets 
in triple systems and on the metallicity difference between 
components with and without planets. In Sect. |S]we sum- 
marize our conclusions. 



2. Observations and data reduction 

2.1. High-resolution spectroscopy with SARG 

SARG, the high-resolution spectrograph of TNG (Gratton 
et al. I200ip . was used to gather high-resolution spec- 
tra of the components of HD132563. The spectrograph is 



equipped with an iodine cell to obtain high-precision radial 
velocities. The spectra were taken with the yellow grism 
(spectral coverage between 4620 and 7920 A) , which is ideal 
for including the full range over which the lines of the io- 
dine cell are present and which at the same time provides 
the highest instrument efficiency. We used The 0.25 arcsec 
slit, which yields a 2-pixel spectral resolution of 144000. 
All the spectra except one (used as template for the RV de- 
termination and for stellar characterization purposes) were 
obtained with the iodine cell inserted into the optical path. 
Integration time was fixed at 900 or 1200 s, to keep er- 
rors due to the lack of knowledge of flux mid-time of each 
exposure smaller than photon noise. After IID132563B was 
selected as a high-priority planet candidate, two consecu- 
tive spectra were acquired to make RV errors smaller. In 
these case, nightly averages were computed and used in 
the following analysis. The templates were obtained with 
an integration time of 3600 s in optimal seeing conditions 
(0.75 arcsec). Median seeing of star-|-iodine spectra, mea- 
sured as FWHM of the spectra across dispersion, is 1.23 
arcsec. While the contamination of the spectra is not of 
particular concern given the projected separation of this 
pair (4.1 arcsec), we followed our standard observing pro- 
cedure that foresees a slit oriented perpendicularly to the 
binary projected separation to minimize the contamination 
of the spectra by the other component. Overall, 53 and 47 
spectra with the iodine cell of IID132563 A and B respec- 
tively were acquired from June 2001 to April 2011. Data 
reduction was performed in a standard way using IRAF. 
Radial velocities were obtained using the AUSTRAL code 
(Endl et al. 1^0001) . as described in Desidera et al. (^0051) . 

2.2. Imaging with AdoptQTNG 

HD 132563 was observed with AdOpt @TNG , the adaptive 
optics module of TNG (Cecconi et al. 2006), with the goal 
of detecting or constraining the mass of the companion re- 
sponsible of the RV variations of HD132563A and to better 
characterize the circumstellar environment of both compo- 
nents. HD 132563 was observed on 20 Jul 2007 (a series of 
22 images) and on 10 Aug 2008 using the Br7 intermediate- 
band filter (three series of 27, 23, and 29 images at differ- 
ent rotation angles to better remove instrumental artifacts) . 
Individual integration times were chosen to avoid detector 
saturation from the bright components, as we are aiming to 
look for companions at small angular separation, where the 
responsible for the RV variations of HD132563a is expected 
(see Sect. 15. 2p . Pointing jitter was performed for sky sub- 
traction purposes. The primary was used as reference star 
for adaptive optics. 

The instrument has a field of view of about 44 x 44 
arcsec, with a pixel scale of 0.0437 "/pixel. Plate scale and 
absolute detector orientation were derived by us when ob- 
serving the systems with long term trends and planet can- 
didates of our survey (see Desidera et al. 120101 for some 
preliminary results). 

Data reduction was performed by first correcting for de- 
tector crosstalk using routines developed for this purpose at 
TNGQ and then performing standard image preprocessing 
(fiat-fielding, bad pixels, and sky corrections) in the IRAF 
environment. 



www . tng . iac . es/instruments/nics/f iles/crt_nics7 . f 
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Fig. 1. Image of HD132563 system obtained with 
AdOpt@TNG. Left and central panels: HD 132563 A and 
B, respectively. Right panel: image obtained subtracting 
the companion PSF for each component. The red circles 
show the 1 X/ D limit for each component. 



be small, thanks to the similar colors of the stars. We con- 
servatively adopt an error of 0.010 mag. 

Table 1. Relative astrometry of the components of HD 
132563 from the observations with AdOpt@TNG. 



Epoch 


e 

deg 


P 
arcsec 


2007.548 
2008.608 


276.69±0.12 
276.95±0.16 


4.109±0.015 
4.130±0.015 



The two components are within the isoplanatic angle 
and then characterized by a very similar PSF (Figure [!}. 
This allows us to achieve much better detection limits on 
difference images, subtracting the PSF of one component 
from the other. 

Images with the same rotation angle were reduced to- 
gether to obtain a single image for each different rotation. 
As a first step we created two separate images by masking 
one of the components. On these images we found the po- 
sition of the center of the component that is not masked 
by applying a two-dimensional Gaussian fit procedure. We 
then subtracted the stellar profiles of both the components 
by subtracting the standard deviation in circular annuli 
around the position of the star. We then applied a proce- 
dure to remove bad and hot pixels that simply substitute 
the value of a pixel with the median of the eight surround- 
ing ones if its value is higher than the value of the median 
plus four times the standard deviation on these pixels. All 
the images were then shifted so that the position of the 
stars is always the same. Finally, we applied a high-pass 
filtering procedure that substitutes the value of each pixel 
with the median of a square of n x n pixels around it, where 
however, we excluded the central pixels to avoid subtracting 
eventual objects. The value of n can be changed according 
to the distance from the central star. At the end of the pro- 
cedure, we made a mean of the the two masked images to 
get the final reduced image. We then made a median on the 
datacubes containing all the images with the same rotation. 
We did discard the images with a value of the FWHM that 
was higher than the value of the median plus the standard 
deviation of the distribution. When images taken at differ- 
ent rotation angles are available, the resulting images from 
this procedure were then rotated and summed together to 
enhance the background objects in the FOV. To further re- 
duce the speckle noise in the final image we subtracted each 
other the square regions composed by 51 x 51 pixels around 
the position of the two components. Before the subtraction 
the two subimages were normalized, taking their total fiux 
into account. Figure [T] shows the image after subtracting of 
the two components of the HD132563 system. Table [T] lists 
our measurements of the relative astrometry between the 
components. 

The magnitude difference in the Br-^ band was obtained 
through aperture photometry. Internal errors are within 
0.005 mag. Additional errors due to transformation on stan- 
dard system cannot be determined from our observations. 
However, we expect the effect on magnitude difference to 



3. Stellar properties 

HD 132563 (HIP 73261, ADS 9461) is a wide pair of main 
sequence, solar-type stars. The physical associations of the 
components was confirmed by the Hipparcos and histori- 
cal astrometry, the small RV difference between the com- 
ponents, and the spectral properties of the objects (see be- 
low). The main properties of the components of HD 132563 
are summarized in Table [5J An additional component (ADS 
9461 C) at about 1 arcmin is listed in several binary star 
catalogs (e.g. CCDM, WDS). However it is not physically 
associated to the system (different proper motion with re- 
spect to HD 132563AB; Roeser et al. i2008)) . We note that 
the photometry of the individual components likely has 
larger errors than those quoted in TableJU because of blend- 
ing effects. 

3.1. Chemical composition 

The iron abundance of the components of HD132563 was 
determined in Desidera et al. (j2004p . The stars result 
moderately metal-poor ([Fe/H]=-0.19). There is no signif- 
icant iron abundance difference between the components 
(A[Fe/H]= 0.013 ± 0.012). The implication of this rcsuft is 
discussed further in Sect. 17. 2[ The abundance of a elements 
compared to iron is close to solar, as expected for thin-disk 
stars. 

Lithium abundance was derived through spectral syn- 
thesis of the 6708 A Li doublet, adopting the atmospheric 
parameters from Desidera et al. (|2004p . We found log N^i = 
2.7 ± 0.1 for HD 132563A and log NL^ = 2.7 ± 0.1 for HD 
132563B. 



3.2. Isochrone fitting and stellar masses 

After we adopted the stellar parameters from the abun- 
dance analysis, the stellar models by Girardi et al. (200^ 
using the web interface pararrO: Da Silva et al. I2006[) 
yield stehar masses of 1.071 ± 0.045 for HD 132563A and 
0.991 ± 0.046 Mq for HD 132563B and ages of 3.6 ± 3.0 
and 4.4 ± 3.9 Gyr. However, these estimates do not con- 
sider the much smaller error bars on the difference of mag- 
nitude, temperature, and metallicity between the compo- 
nents. The magnitude difference A.V, coupled with the ef- 
fective temperature difference ATog, indicates that the pri- 
mary is slightly evolved. After adopting the spectroscopic 
effective temperatures and metallicities, the age is tightly 



|http : //stev . oapd . inaf . it/paraiii| 
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Table 2. Stellar properties of the components of HD 
132563; where one value is listed at the center of the table 
means the measurement refers to joined A+B components 
or to the difference between the components. 



Parameter 



HD 132563 A HD 132563 B Ref. 



a (2000) 


14 58 21.519 14 58 21.136 


1 


S (2000) 


+44 02 35.33 +44 02 35.87 


1 


fia (mas/yr) 


-59.86 ± 1.41 


2 


^s (mas/yr) 


-70.12 ± 1.19 


2 


RV (km/s) 


-6.68+ 0.20=" -6.05+0.20 


3 


RV (km/s) 


-6.5+1.0*^ 


3 


TV (mas) 


10.41 + 1.45 


2 


d(pc) 


96 + 13 


2 


U (km/s) 


8.8+1.1 


3 


V (km/s) 


-39.3+4.9 


3 


W (km/s) 


13.5+2.8 


3 


rmin (kpc) 


6.30+0.31 


3 


rmax (kpc) 


8.594+0.024 


3 


Zmax (kpc) 


0.277+0.039 


3 


e 


0.154+0.023 


3 


V 


8.966+0.007 9.472+0.010 


4 


AV 


0.505+0.018 


4 


B-V 


0.560+0.023 


1 


B-V 


0.54+.03 0.57+0.05 


4 


V-I 


0.63+0.02 


1 


Hp scatter 


0.016 


1 


<V2A/ass 


7.731+0.054 


5 


H2Masa 


7.576+0.033 


5 


K2Mass 


7.431+0.21 7.494+0.049 


5 


AK 


0.467+0.010 


3 


NUV magnitude 


13.465+0.003 


6 


FUV magnitude 


19.27+0.08 


6 


Mv 


4.05+0.30 4.56+0.30 


3 


Teff (K) 


6168+100'= 5985+100 


4 


AT^ffiA-B) (K) 


184+12 


4 


logff 


4.15^= 4.27 


4 


V sin i (km/s) 


2.0+0.8" 1.5+0.5 


3 


[Fe/H] 


-0.18 + 0.10== -0.19 + 0.10 


4 


A[Fe/R](A~B) 


0.012 + 0.013 


4 


EW Li 6707 A 


30.2'= 42.9 


3 


log NLi 


2.r= 2.7 


3 


Mass(M0) 


1.081+ 0.010^= 1.010 + 0.010 


3 


Age (Gyr) 


~5 


3 



References: 1 Hipparcos (ESA 19,97J; 2 van Leeuwen (2007); 3 
Th is Pa per; 4 Desidera et al. (i2004|l 5 2MASS (Skrutskie et 
al. lMI^ : 6 GALEX (Martin et al. lM?g|) . 

"" Instantaneous RV of HD 132563 Aa 

^ Approximate system RV taking into account all the known 

components in the system 
'^ The value refers to HD132563Aa component only. 




constrained at 5.5 Gyr (Fig. [2]). Uncertainties in tempera- 
ture scale and absolute abundance implies an error of about 
1.5 Gyr on stellar age. The stellar masses determined from 



Fig. 2. Isochrone fitting of HD132563 system. The Girardi 
et al. ([2002)1 3.0, 5.0, 5.5, and 6.0 Gyr Z=0.013 (the nominal 
metallicity of the system) isochrones are shown, together 
with the observational data. The larger error bars (dotted 
lines) refer to the absolute error (dominated by error on 
distance); the smaller error bars (continuous lines) refer to 
errors on the difference between the components. 



the best-fit isochrontS are 1.081 and 1.010 M© for HD 
132563 A and B, respectively. 



3.3. Activity and rotation 

Projected rotational velocities were measured by fast 
Fourier transform (FFT) of absorption profiles, which were 
computed by the cross correlation function (CCF) of the 
stellar spectra involved in the line bisector study. About 
30-35 unblended lines, mainly Fe in the wavelength range 
5230-6050 A, were selected when building the CCF. The 
CCF profiles (i.e. the observed profiles from the stars) were 
made symmetric by mirroring one of its halves, aiming to 
reduce the noise of the FFT. A new profile was calculated 
by the convolution of a macroturbulence profile (Gaussian) 
and a rotational one, to compare the FFTs of the symmetric 
and the calculated (model) profiles. Macroturbulence was 
estimated from effective temperature using the standard re- 
lation by Valenti & Fischer (|2005p . The projected rotation 
velocity v sin i corresponds to the case where the first mini- 
mum of the FFT from the calculated profile coincides with 
the first minimum of the FFT from the symmetric one, as 
described by Gray (|1992p . The projected rotation veloci- 
ties of HD 132563 A and B result vsini = 2.0 ± 0.8 and 
1.5 ± 0.5 km/s, respectively. 

A direct measurement of magnetic activity level would 
be very useful for interpreting the RV variations in terms of 
Keplerian motion or activity jitter. Unfortunately, it is not 
possible to simultaneously observe the region that us rich in 
iodine lines and the Call H&K chromospheric lines using 
SARG. An indirect estimate in the chromospheric emis- 
sion can be derived from the measurement of UV excess 
in GALEX photometry. Adopting the joined B and V pho- 
tometry from Hipparcos and the FUV GALEX magnitude 
of 19.27 ± 0.08, we get an UV excess of -0.94 mag over the 
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Fig. 3. The FFT of HD132563 A and B (upper and lower 
panels) for three absorption profiles: the solid line for the 
observed profiles, the dashed line for the symmetric profiles 
where the FFT appears less noisy, and the dotted line for 
the calculated (model) profiles. 



Table 3. Summary of age estimates of HD132563 compo- 
nents. 



Method 



Age 
Gyr 



Isochrone fitting 5.5±1.5 
Chromospheric emission ~ 1.2 
X-ray emission > 1.3 

Lithium 1-4 

Rotation 2-6 

Kinematic 3-6 



Adopted 



3.4. Absolute radial velocity and Kinematics 

The absolute radial velocity of the components obtained 
through cross-correlation with templates of stars from the 
list of Nidever et al. ([2002) results of -6.68 ±0.20 km/s for 
HD132563 A and -6.05 ±0.2 km/s for B respectively. The 
systemic velocity to be used for the kinematic analysis can- 
not be determined with accuracy because of the unknown 
mass and orbit of the spectroscopic companion around 
HD132563A. From plausible values derived in Sect. 15.21 
we adopt a RV of —6.5 ± 1.0 km/s. The space velocities, 
taking trigonometric parallax and proper motion from van 
Leeuwen ((2007|) into account, result in U,V,W=+8.8± 1.1, 
-39.3 ± 4.9, and 13.5 ± 2.8 km/s. The UVW velocities 
are well outside the kinematic space populated by young 
stars (Monies et al. I200ip . The galactic orbit was derived 
following the procedure described in Barbieri & Gratton 
(|2002L which assumes the Galactic gravitational poten- 
tial by Allen & Santillan I199ip . The numerical integra- 
tion of the galactic motion yields the following parame- 
ters: r™„ = 6.30 ± 0.31 kpc, r^ax = 8.594 ± 0.024 kpc, 
z„^^ = 0.277 ±0.039 kpc, and e = 0.154 ±0.023. The kine- 
matic properties then indicate an age similar to the Sun 
and possibly even slightly older. 



photospheric value, when using the calibration by Smith 
& Redenbaugh (j2010p . Such an excess is what is expected 
for a star with the color of HD132563 and with log Rhk 
about -4.^ (see Fig. 4 in Smith & Redenbaugh [2{?TU|) . 
The corresponding age using the calibration of Mamajek 
& Hillenbrand (2008) is 1.2 Gyr, with a large uncertainty. 
The X-ray emission of the object was not detected by 
ROSAT (Voges et aL [2000|) . The upper hmit of logLx < 
28.6 (assuming equal X-ray luminosity for the two com- 
ponents and neglecting the tertiary component around the 
primary, see below) corresponds to an age older than about 
1.3 Gyr using the calibration of Mamajek & Hillenbrand 
(|2008p . The star also does not have large photometric 
variability (0.016 mag dispersion from Hipparcos measure- 
ments, which refer to the joined A+B components and is 
flagged as "duplicity induced variability", so that the in- 
trinsic variability is likely below this value). 



* An alternative explanation for the small UV excess is rep- 
resented by a hot white dwarf in the system, possibly the spec- 
troscopic companion of HD132563A. However, there are several 
reasons, discussed in Sect. [5l to consider more likely that the 
spectroscopic companion is an M or K dwarf rather than a white 
dwarf. 



3.5. Stellar age 

The isochrone age of 5.5 ± 1.5 Gyr is fully consistent with 
the kinematic properties of the system, the slow projected 
rotational velocity, and the X-ray non-detection. Lithium 
and chromospheric activity instead favor a younger age of 
about 1-3 Gyr. Considering the lower error in the isochrone 
fitting, we adopt 5 Gyr as the most probable stellar age 
(Table O. 



4. A giant planet around HD132563B 

4.1. Radial velocity variations of HD132563B 

The RVs of HD132563B (Table El show significant varia- 
tions (r.m.s. 23.4 m/s) compared to internal errors (10.3 
m/s). Internal errors are larger than those for the typical 
targets of the SARG planet search because the object is at 
the faint end of our sample. The Lomb-Scargle periodogram 
fScargle ri982p shows a well-defined peak at about 1500 d. 
Bootstrap simulations performed by scrambling the RVs 
yield a false alarm probability of 5/10000 (Fig. H]). This 
is low enough to provide confidence that the detected RV 
periodicity is real. 
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Fig. 4. Upper panel: radial velocities of HD 132563 B, over- 
plotted is the Keplerian best fit. Central panel: radial ve- 
locities phased to the best-fit period. Lower panel: Lomb- 
Scargle periodogram of the radial velocities. Overplottcd 
are the various confidence levels as derived from the boot- 
strap analysis. 



4.2. Origin of RV variations 

Before inferring the Keplerian nature of the observed vari- 
ations we have to exclude alternative possibilities such as 
magnetic activity. An activity jitter of 20 m/s is usually ob- 
served in stars of similar color with log R'^k ^ —4.4-; — 4.5 
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Fig. 5. Upper panel: bisector velocity span (BVS) vs radial 
velocity for the spectra of HD132563B, showing no corre- 
lation. Lower panel: Line bisectors of individual spectra of 
HD132563B. The parts of the line profile adopted for the 
derivation of BVS are shown as horizontal lines. 



and age of about 400 Myr. As discussed in Sect.|31 we do not 
have any indication of such a strong magnetic activity of 
the star or such a young age. Furthermore the observed RV 
period is much longer than the expected rotational period 
of a solar type star. 

Finally, activity-induced variations are usually accom- 
panied by changes in line profiles correlated with those of 
radial velocities, even if this diagnostic is less senstive for 
slowly rotating stars such as HD132563B; see Santos et 
al. (pbn5|) and Desort et al. POITTD . To this aim, we per- 
formed the analysis of line bisectors using the technique 
described in Martinez Fiorenzano et al. (|2005p . As shown 
in Fig [SI there is no significant correlation between RV and 
bisector velocity span (BVS). 

Another potential source of spurious RV variations in 
multiple systems is represented by contamination of the 
spectra by the light of the other components. This is not 
of concern in the present case considering the projected 
separation of about 4.1 arcsec. A few spectra taken with a 
seeing of about 2 arcsec show no indication of significant 
contamination neither as outliers in the RV curve nor in 
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Table 4. Orbital parameters and results of fitting for RV 
of HD132563B. 



Parameter 




Value 


Period (d) 




1544±34 


K (m/s) 




26.7±2.2 


e 




0.22±0.09 


uj (deg) 




158±35 


TO 




2593±148 


msini (Afj) 




1.49±0.09 


a (AU) 




2.62±0.04 


r.m.s. res (m 


/s) 


12.7 



the line profile and there is no correlation between RV and 
seeing. 

We then conclude that the most likely explanation for 
the observed RV variation is the presence of a low-mass 
object, HD132563Bb, whose projected mass is clearly in 
the planetary regime. 

4.3. Orbital parameters 

Orbital parameters to characterize the detected RV vari- 
ations of HD132563B were determined by using an our 
own code, which is based on a Levenberg-Marquardt least- 
squares fit of the RVs. An additional error of 7 m/s (repre- 
senting activity jitter and possible unrecognized systematic 
errors) was added in quadrature to internal errors when per- 
forming the fit. The best-fit parameters are listed in Table 
0] and the best-fit solution plotted in Fig. U) The errors in 
the orbital parameters were derived by simulating synthetic 
datasets that take the adopted jitter and performing orbital 
fitting into account for each fake RV series. The r.m.s. of 
residuals from the fit are 12.7 m/s, slightly larger than the 
internal errors. A circular orbit fitting yields slighly larger 
residuals (13.5 m/s), indicating that the value eccentricity 
is rather uncertain. From the stellar mass 1.010 Mq, we 
derive a projected mass of 1.49 Mj and a semimajor axis 
of 2.62 AU for HD132563Bb. 

Analysis of residuals from the best-fit orbit does not 
reveal any significant periodicities (Fig|6]). There is a pos- 
sible long-period modulation that is not highly significant 
but that calls for the continuation of the RV monitoring to 
look whether it is real and, in this case, if it might be due 
to an additional planet in outer orbit. The highest peak in 
the periodogram of residuals in the range 10-30 days is at 
16.4 days, which is different from the second highest peak 
in periodogram of the originl RVs (15.3 days). Such a pe- 
riodicity might be related to stellar rotation but can just 
come from noise fluctuations. 



5. A stellar companion in highly eccentric orbit 
around HD132563A 

5.1. Radial velocity variations of HD132563A 

The RV of HD132563A showed a nearly linear trend for a 
few years. The RV maximum was reached in 2006. After 
this, an RV decrease with increasingly steep slope was ob- 
served (Fig. [7]). The peak-to- valley RV variations are now 
more than 1 km/s, placing the companion to HD132563A 




Fig. 6. Upper panel: Residuals from best fit orbit Lower 
panel: Lomb-Scargle periodogram of residuals from formal 
Keplerian best fit of the radial velocities of HD132563B. 



well outside the planetary mass regimqj. The orbit is highly 
eccentric. Residuals from the formal best-fit orbit do not 
show any significant periodicities. The dispersion of the 
residuals is 14.1 m/s, which are significantly larger than 
the internal errors (8.5 m/s), leaving room for a jitter of 
about 11 m/fl 

Because only a portion of the orbit of HD 132563 A 
was observed, a wide variety of orbital solutions and com- 
panion masses are compatible with the RV data (Sect. 15. 2p 
However, additional, available data allow us to further con- 
strain the possible parameters of the companion. They are 
described below individually, and conclusions on the possi- 
ble mass and orbit are presented in Sect. 15.81 



'' Following standard notation, we identify the primary of the 
spectroscopic pair as HD 132563 Aa, contributing to a large frac- 
tion of the observed flux, and its spectroscopic companion as 
HD132563Ab. 

^ An additional error in the RVs might also be expected if 
the contribution of HD132563Ab to the integrated spectrum is 
not negligible. In this case, the velocity difference between the 
spectroscopic components changes with time, and the template, 
which is unique for all the observations, might not be fully rep- 
resentative of the true stellar spectrum at the various epochs, 
causing some degradation in the modeling of the composite io- 
dine-|-star spectrum. 
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Table 5. Possible orbital parameters of HD132563A, con- 
sidering the orbital solutions that give r.m.s of residuals 
< 15.7 m/s. 
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Fig. 7. Radial velocities of HD 132563A. Overplotted the 
Keplerian best fit. 



5.2. Montecarlo simulation 

To take the various constraints on the mass and orbital pa- 
rameters of the spectroscopic companion into account we 
performed a dedicated MonteCarlo simulation. The pro- 
cedure first generates a set of companion masses and all 
the orbital parameters. We established an upper limit on 
the semimajor axis at 40 AU (corresponding to a period of 
about 200 yr). Wider orbits appear unlileky because typ- 
ically the inner pair in triple systems has an orbit that 
is much closer than the outer component, well within the 
dynamical stability region and the probability of catching 
the phase of steep RV variation at periastron passage in 
an RV monitoring of about ten years is lower than 5%. 
Furthermore, the relative astrometry does not indicate any 
periodicities that are a few times the time baseline of the 
observations, which would be expected because the astro- 
metric amplitude is larger at longer periods for given com- 
panion mass, and best-fit solutions in our simulation are 
confined to shorter periods (see below) . We also fix the up- 
per limit of the companion mass at 0.8 Mq, from the lack of 
signatures in our spectra fSect. 1^751) . From each generated 
input set, we derived the corresponding RV variations at the 
time of SARG observations, the projected separation and 
luminosity contrast at the time of AdOpt observations, and 
the astrometric amplitude for comparison with the obser- 
vational data. This allowed us to identify the combination 
of parameters that produces an RV curve compatible with 
the observed one (r.m.s. of residuals < 18.5 m/s, 2cr limit 
as compared to the residuals from the formal best-fit orbit 
using F test; Fig.jT]). Figure |S] shows the results of the sim- 
ulations: there is a broad range of orbital parameters that 
are compatible with the RV signature, which becomes nar- 
rower when considering only the best-fit solutions {la limit; 
r.m.s. of residuals < 15.7 m/s). The orbit results highly ec- 
centric (e > 0.65), with periastron passage expected within 
the end of 2012. Table [5] lists the most probable orbital 
parameters as resulting from our simulations. 

5.3. Spectral properties and line bisectors 

We analyzed the spectra of HD 132563 A and the result- 
ing line profiles looking for signatures of the spectroscopic 



Parameter 


Mean 


Median 


r.m.s. 


a(AU) 


14.8 


12.9 


6.3 


P(yr) 


47 


37 


30 


e 


0.860 


0.874 


0.065 


K (m/s) 


5460 


4887 


2995 


msini(AfQ) 


0.40 


0.38 


0.18 


mass (Mq) 


0.49 


0.51 


0.20 


p (arcsec) 


0.145 


0.146 


0.018 


uj (deg) 


160.2 


160.0 


5.9 


TO (yr) 


2012.34 


2012.34 


0.20 


7 (m/s) 


-1625 


-1613 


288 


astrom. sign, (mas) 


51 


40 


43 



companion on the spectra (Fig. ^. There is no significant 
correlation between RV and BVS and the individual line 
profiles typically appear clean without indicating the spec- 
tral signatures of the spectroscopic companion. As a rough 
guess, we estimate that secondaries contributing more than 
10% of the light in the T^band are excluded by the line pro- 
file analysis. The absolute shape of the line bisector, less 'C 
shaped' than HD132563B, might indicate of a small con- 
tamination of the line profile at negative RV (where we 
expect the contribution of the spectroscopic companion), 
at the level of a few %. 

5.4. Magnitude difference between the components 



As shown in Sect. 13.21 a fairly good isochrone fitting is 
obtained by exploiting the small errors on the V magnitude 
and temperature differences. The K mag difference, 0.467 ± 
0.010, is not compatible with what is expected from the best 
fit isochrone. The magnitude difference is larger by about 
0.09 mag. 

Because we know that there is an additional compo- 
nent in the system, HD132563Aa, we derive its properties 
assuming that it is the cause of the K band excess and 
that the other two components follow the best-fit isochrone. 
We found Mk=2.76, 3.14, and 5.45 for HD132563Aa, 
HD132563B, and HD132563Ab, respectiv ely. U sing the 
mass- luminosity relation by Delfosse et al. (|2000p . the cor- 
responding mass of the spectroscopic companion would be 
0.56 A/©. The relative contribution of HD132563Ab in V 
band would be smaller, about 0.7%, with little effect on 
the isochrone fitting. 

5.5. Direct imaging 

No companions close to either HD132563 A or B were de- 
tected in our images taken with AdOpt@TNG (Sect. [^ . 
The procedure of subtracting the PSF of the other compo- 
nent leaves residuals that cannot be eliminated by optimiz- 
ing the position shift between the components. This indi- 
cates a slightly different PSF between the two components. 
A companion contributing less than 10% at a very small 
projected separation might produce such feature, which 
might also, however, be related to instrumental effects. 

To estimate detection limits, we considered the disper- 
sion of fluxes of pixels in circular annuli around the cal- 
culated position of the star, with steps of one pixel. The 
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Fig. 8. Results of the MonteCarlo simulation, showing the possible parameters that are compatible with the observed 
RV curve of HD132563A. Starting from top left: eccentricity, RV semiamplitude, companion mass, time of periastron 
passage, longitude of periastron, RV difference between the components vs semimajor axis; RV difference between the 
components vs companion mass; companion mass vs projected separation (at epoch 2007.548); astrometric amplitude vs 
orbital period. Blue circles: solution with r.m.s. of residuals < 15.7 m/s; Green circles (light gray in BW): solution with 
r.m.s. of residuals 15.7 — 18.5 m/s. The red solid line in the mass vs RV difference panel is the limit that allows the wide 
pair to be bound. The red solid line in the companion mass vs projected separation panel shows AdOpt@TNG detection 
limits. 



detection limits were set at five times the standard devia- 
tion in each annulus. During these calculations the pixels at 
less than 30 pixels from the other component were masked. 
In Fig. [TU] we display the contrast obtained around the A 
component and the corresponding mass detection limits. 
To translate the contrast limit on companion masses, we 
adopt the K magnitude of the components from 2MASS 



(Skrutskie et al. '2006111 and use the mass-luminosity rela- 
tion by Delfosse et al. (2000) for companions with Mk < 9 
and Chabrier et al. (|20D0| models at fainter magnitudes 



(using the adopted stellar age). 

The expected projected separation at the epoch of our 
observations is between 0.09 and 0.18 arcsec (Fig. [8]) at 



^ 2MASS photometry for close pairs of bright stars is rather 
poor, with the error exceeding the nominal ones because of 
blending. 
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Fig. 9. Upper panel: Bisector velocity span (BVS) vs radial 
velocity for the spectra of HD132563A. No correlation is 
present. Lower panel: Line bisectors of individual spectra 
of HD132563A. The parts of the line profile adopted for 
deriving BVS are shown as horizontal lines. 
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Fig. 10. K band detection limits for observations taken 
with AdOpt@TNG in July 2007 (dashed hne) and August 
2008 (continuous line). A few individual values of masses 
corresponding to contrast values are also shown. 



2007.548 and smaller by about 10-15% one year later, as 
the pair is approaching periastron passage. The contrast 
limit derived at such separations corresponds to about 
0.6 Mi^ Detection limits at separations larger than 1.5 
arcsec (background-limited region) is close to the stel- 
lar/substellar boundary. 

The mentioned mass limits apply to main sequence 
stars. White dwarfs are faint in K band (Bergeron et 
al. I200ip and are not excluded by our direct imaging 
data. However, a white dwarf appears unlikely as close 
companion. HD 132563A is a close enough binary that 
the (current) primary HD132563Aa should have accreted 
a significant amount of material with the signature of 
AGB nucleosyhthesis, becoming a barium dwarf, if its 
companion HD132563Ab has evolved through the AGB 
phase and is now a white dwarf (see McClure et al. 119801 
Frankowski & Jorissen [2007J) . However, the barium abun- 
dance of HD132563A is identical within 0.01 dex to that 
of HD132563B. The small FUV excess of the system can 
be explained by a moderate chromospheric activity of the 
stars. 



5.6. Relative astrometry 

The relative astrometry between the components of 
HD132563 can be used to further constrain the proper- 
ties of the companion around HD 132563A. Available as- 
trometric measurements in the Washington Double Star 
Catalog (Mason et al. [^UOTl version 4 Apr. 2011, kindly 
provided to us by B. Mason) are shown in Fig. [TTJ The as- 
trometric time series over a span of about 180 years clearly 
show the signature of the orbital motion. As no curva- 
ture is detected, we assume this is due to the motion of 
the wide pair and not to the spectroscopic companion to 
HD132563A (otherwise the period of the companion would 
be very long, making the probability of catching the ob- 
ject at periastron passage very low). The linear best fit 
yields slopes of dp/dt — —0.00353 ± 0.00036 arcsec/yr and 
dO/dt = -0.0185 ± 0.0033 deg/yr. 

Figure [TH shows the relative astrometry of HD132563 
in the last 20 years considering only high-quality mea- 
surements: Hipparcos, speckle interferometry (Douglass et 
al.Hnnni Scardla et al. 2007; Losse l^IITUI and our own AdOpt 
measurements from Table [1]). While the details of the resid- 
uals pattern depend on the adopted long-term slope, which 
has some uncertainty, the measurement by Douglass et 
al. (2000) is well outside of the linear fit obtained consider- 
ing all epochs back to 1830. It is plausible that we are seeing 
the astrometric impact of HD132563Ab on HD132563Aa. If 
this is the case, a lower limit to the amplitude of the astro- 
metric motion would be about 50 mas. This is compatible 
with the astrometric signature of a 0.56 AIq companion, 
as suggested by the K band excess, and the most probable 
orbital elements from RV. 



5. 7. Radial velocity difference between the components 

As HD132563A is approaching the periastron of the spec- 
troscopic orbit, its RV becomes more negative, and the RV 



Detection limits are uncertain at small projected sepa- 
rations, where there are few points in the circular annuli. 
Furthermore, a moderately bright companion can affect the cen- 
tering procedure. 
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Fig. 11. Relative motion of HD 132563. Continuous line is 
the linear fitting to the data. 
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Fig. 12. Recent high-quality relative astrometry of 
HD132563 system. Continuous line is the linear fitting 
to the data shown in Fig. 1111 The measurements from 
Douglass et al. (,2000|1 and AdOpt@TNG were averaged. 



difference between the components increases. However, for 
increasing values of the RV difference between HD 132563 
A and B (because of a high value of RV amplitude of the 
spectroscopic variations) it would become no longer pos- 
sible to ensure that the two wide components are bound. 
We derived such a limit using the approach presented by 
Hauser & Marcy PM^I) . 

We first obtained the RV difference between the com- 
ponents of HD132563 by measuring the RV of HD132563B 
with the stellar template of HD132563A. We found an RV 



difference between the mean of the RVs of each compo- 
nent of ARV{A - B) ^ -747 ± 35 m/s. The mean RV 
of HD132563B is very close to the center of mass of the 
star-|-planet system, while the center of mass of HD132563A 
is unknown because of the still undetermined spectroscopic 
orbit. However, from the orbital solutions compatible with 
the RV data we can infer that the absolute RV difference 
between the components should be between —900 to —2100 
m/s, with a dependence on the mass of the spectroscopic 
companion (and semimajor axis, see Fig. [5] Sect. 15.2]) . 

From the relative astrometry of the components, assum- 
ing that what we observe is their relative orbital motion, 
we derived the relative positions and velocities between the 
components on the plane of the sky {x,y,Vx,Vy)- These 
result a; = 52 AU, y = -397 AU, v^ = -814 m/s, and 
Vy = 1515 m/s. The RV difference between the components 
yields the velocity along the line of sight (w^). The relative 
position along the line of sight (z) remains unknown. 

The maximum separation for having a bound pair is 



= y/x^ +y^ + zl 



2.0 G iriA TUB 



(1) 



where itia and m^ are the masses of HD132563 A and B, 
respectively, m^ed is the reduced mass, and v the orbital 

velocity v — Jv'^ + Vy + v^ 

The RV differences between A and B that do not al- 
low the system to be bound would be lower than —2.67 
and -3.07 km/s for masses of HD132563Ab of 0.20 and 
0.70 Mq respectively. The actual values of the plausible 
orbital solutions are well above this limit (Fig. [S]). 

5.8. Conclusion on the characteristics of the spectroscopic 
companion 

The results described above does not allow us to fully and 
conclusively determine the mass of HD132563Ab but are 
all consistent with a 0.5 — 0.6 M© companion. Such a mass 
is suggested by the small K band excess and the tenta- 
tive evidence of astrometric signature. It is close to the 
most probable value suggested by the Montecarlo simula- 
tion. Contribution to the integrated flux of HD132563A is 
expected to be at a < 1% level in the visible, with negli- 
gible impact on the various results presented in the paper. 
The RV variations strongly constrain the orbital period, 
eccentricity, longitude of periastron and epoch of perias- 
tron passage. The continuation of the RV monitoring will 
allow a firm determination of these parameters and of the 
projected mass. 

6. Binary orbit 

The approach presented in Sect. 15.71 can also be used to 
put constraints on the orbit of the wide components (see 
Hauser & Marcv [T999p . We use the values of a;, y, v^ and 
Vy given above and explore the orbital parameters result- 
ing from the allowed values of z, the separation between the 
components along the line of sight. Because a range of val- 
ues for the RV difference and the total mass HD 132563A 
are possible, we perform the calculation for three repre- 
sentative pairs of these quantities (taking the correlation 
into account that results from our MonteCarlo simulations, 
see Fig. [5]): |Ai?y|=1.2, 1.4, and 1.6 km/s and mass of 
HD132563Ab of 0.20, 0.45, and 0.70 Mq, respectively. The 
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Fig. 13. Seminiajor axis, eccentricity, periastron, and incli- 
nation of the binary orbit as a function of the separation 
along the line of sight z for various RV differences between 
the components and masses of HD132563Ab. Continuous 
line: RV difference -1.4 km/s and mass 0.45 M©. Dotted 
line: RV difference -1.2 km/s and mass 0.20 Mq. Dashed 
line: RV difference -1.6 km/s and mass 0.70 Mq. 



results are shown in Fig. [T31 while Fig. [T3] shows the cor- 
responding critical semimajor axis for dynamical stability 
{acrit ; Holman & Weigert I1999[FI . The plots of the orbital 
parameters for the three adopted configurations look qual- 
itatively similar. Only highly eccentric orbits are allowed 
for z < 0, with rather close separations at periastron and 
limited regions for dynamical stability around both compo- 
nents. Such orbits of the wide binary are then excluded by 
the presence of the companion orbiting HD132563A, which 
achieves a separation at apoastron of at least 15 AU. The 
most plausible zone is for z between 100 and 700 AU, cor- 
responding to orbits with a - 400 - 3000 AU, e ~ 0.2 - 0.8, 
periastron larger than ^ 100 AU, and acrit larger than 30 
AU for both components of the wide pair. 



7. Discussion 

7.1. The planet around HD132563B 

The planet around HD132563B is a rather typical giant 
planet, lying in a well populated region in planetary mass 
vs orbital period or semimajor axis diagram. The eccen- 
tricity is moderate and also not unusual for planets at this 
separation. HD 132563B has a lower mctallicity with re- 
spect to the bulk of planet hosts, but its metal deficiency 
cannot be considered really peculiar. 

The orbit of the HD132563 system is still undetermined, 
so we did not perform a detailed derivation of the dynam- 
ical stability of the planet. However, the presence of the 



The equation to derive acrit is defined only for eccentricities 
below 0.8, and therefore the results shown in Fig. [TJ] are not 
accurate for the orbits with extreme eccentricities. 



Fig. 14. Critical semimajor axis for dynamical stability for 
the various wide binary orbits shown in Fig. [T3] (upper 
panel: HD132563B; lower panel: HD132563A). Line styles 
as in Fig. [T3l 



HD132563A pair, with a separation several times larger 
than the planet from its host star, guarantees that the 
planet orbit is stable. 

7.2. Abundance difference between binary components with 
and without planets 

As shown in Sect. [31 there are no significant abundance 
differences between the components of HD 132563 (A 
[Fe/H]= -|-0.012±0.013). Our systematic searches for abun- 
dance differences in binary systems with similar compo- 
nents (Desidera et al. 120041 12006|) allowed us to infer that 
large abundance differences are rare. Only one case with an 
abundance difference larger than 0.10 dex was found (HIP 
64030) but the origin of the abundance variations might 
be unrelated to accretion of planetary material because 
the primary of this system is a blue straggler (Desidera 
et. l2007|) . Also, abundance analysis of stars in open clusters 
did not reveal robust evidence of members with enriched 
abundances (e.g. Paulson et al. 120031 Shen et al. I2005p . As 
we have only one planet known in the sample of Desidera 
et al. 12004! (HD 132563, presented here) and the samples 
of Desidera et al. 12006! and open cluster members were not 
systematically searched for planets, it is possible that the 
lack of abundance anomalies is due to the paucity of planets 
in our sample. We then consider in Table [5] the abundance 
difference in binaries with similar components for which 
giant planets have been detected. While optimized differ- 
ential abundance analysis was only performed for 16 Cyg 
and HD132563, there are no indications of large (> 0.10 
dex) abundance differences in any of the seven systems for 
which data are available. These results indicate that large 
metallicity enhancements comparable to the typical mctal- 
licity difference between stars with and without giant plan- 
ets (e.g. Gonzalez !2006p are rare. We then provide further 
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support to the idea that high metallicity favors the forma- 
tion of giant planets around solar-type stars (e.g. Santos et 
al. 1^004; Fischer & Valenti HOOSl Mordasini et al. [2005)1 . 

Some of the pairs in Table [B] have rather cool compo- 
nents for which the engulfment of fairly large amounts of 
rocky material are required to produce detectable abun- 
dance changes. The case of HD132563 is very favorable 
thanks to the small errors on abundance difference and the 
rather thin convective zone of its components. Following the 
approach by Murray et al. (2001), the masses of the convec- 
tive zones of HD132563 A and B are 0.016 and 0.018 Mq, 
respectively. The 1-sigma error on our abundance differ- 
ence corresponds to the accretion of just 0.13 of iron for A 
and 0.14 M^ for B. These quantities are smaller than the 
amount of rocky material that should have been accreted by 
the Sun during its main sequence lifetime (0.4 Mg of iron 
corresponding to 5 M® of meteoritic material), according 
to Murray et al. ((200T|l . 

Recent results suggest that the Sun has a peculiar abun- 
dance pattern with respect to solar analogs, with a de- 
pletion of refractory elements relative to volatile elements 
(M elende z et al. [MM Gonzalez et al. [20T01 Ramirez et 
al. I2010p . possibly linked to the formation of terrestrial 
planets. A more detailed evaluation of the full relative abun- 
dance pattern of HD132563 components, will be performed 
in a forthcoming paper, by studying elememts with differ- 
ent condensation temperature and exploring these issues. 



7.3. Planets in triple systems 

The planet orbiting HD132563B is one of the few known in 
triple systems. Table [7] lists the current census of planets 
in triple systems. In all these systems the planet orbits the 
isolated component of the hierarchical triple, with a close 
stellar pair at a much larger separation. At first order, the 
dynamical effects of the distant pair can be approximated 
by a single star with the sum of the masses of individual 
components. After deriving the critical semimajor axis for 
dynamical stability (Holman & Wiegert 11999^ in this way, 
there are no systems for which the outer pair is found to 
have a strong impact on the planetary region. Another sim- 
ilar system not included in Table [7] is HD41004, for which 
the companion of the no-planet host is in the substcllar 
regime (msini = 18 Mj; Zucker et al. I2004p . 

A different system architecture might be represented 
by WASP-8=CCDM 23596-3502. This solar type star 
has an M dwarf companion (CCDM 23596-3502B; mass 
about 0.5 Mq) at 400 AU projected separation (Queloz 
et al. I2010p . Besides the transiting planet on a strongly in- 
clined and eccentric orbit, WASP-8 shows a long term drift 
of 58 m/s/yr, without any significant indication of curva- 
ture over two years of observations. Continuation of the 
observation will reveal whether such a companion is a mas- 
sive planet, a brown dwarf, or another stellar object in the 
system. 

For five of the eight hierarchical systems in Table [71 the 
mass of the planet host is lower by more than 10% than the 
sum of the masses of the wide stellar companions. These are 
rather special cases in the current census of exoplanets in 
binaries, for which in nearly all cases the planet host is the 
most massive star in the system. This is likely from a com- 
bination of selection effects. In most cases only the stellar 



primary is under RV monitoring [3, and a significant frac- 
tion of the secondaries are low-mass stars for which there 
is evidence of a lower frequency of giant planets (Endl et 
al. 2006, Johnson et al. I2007p . Solar-type companions of 
early type stars would allow investigating the role of the 
mass ratio on the occurrence of planets. Hierarchical triples 
represent an interesting alternative to such a program, pro- 
viding a sample of stars that are affected by the dynamical 
influence of the equivalent of a more massive companion. 

While a detailed study of the planets in triple systems 
is beyond the scope of this paper, it seems that hierarchi- 
cal triple systems do not represent a hostile environment 
for planet formation around the isolated component, re- 
gardless of the mass ratio between the planet host and the 
sum of the masses of the other components. In the sam- 
ple of our survey, we have detected nine additional com- 
ponents that are likely stellar (19% of the pairs), and the 
first planet was detected in one of such. In the sample of 
planets in binaries studied in Desidera & Barbieri (|2007p . 
there are four triple systems among the 35 considered in 
the study (11%). This is lower than the fraction of triple 
systems with respect to binaries by about a factor of two 
fTokovinin l2008p . but the multiplicity estimate of the com- 
panion of planet host is underestimated (e.g. no dedicated 
RV monitoring looking for short-period companions in sev- 
eral cases). In the sample of Bonavita & Desidera (|20071 
with updates in Bonavita et al. I2010p . there are 25 triple 
systems. For 12 of these systems, the component under RV 
monitoring is a member of the close pair, with no planets 
detected around them. This is consistent with the paucity 
of planets in close binaries (see Bonavita & Desidera 120071 
Eggenberger l2010| ): in these cases, the dynamical effects on 
the target are dominated by close companion. The other 
13 svstemg^ are those for which the isolated component is 
under RV monitoring, with three planets detected (those 
orbiting HD 40979, HD 178911B, and HD 196050), leading 
to a frequency of planet of 23%. Selection biases certainly 
play a role in such a large frequency: dedicated deep imag- 
ing observations such as those that allowed the binarity of 
the companions of HD 40979 and HD 196050 to be identi- 
fied were not performed for the majority of the stars in the 
UD sample. If we consider only the triple systems in the UD 
sample identified independently on the planet detection we 
find a planet frequency of 9% (1/11, HD 178911B), fully 
compatible with what is found for planets around single 
stars and components of wide binaries. This result further 
supports that isolated components in triple systems are not 
devoid of planetary companions. 

In all but one stellar triple with planets, the separa- 
tion of the stellar pair is larger than the planet semimajor 
axis. While selection effects certainly play a role (no spec- 
troscopic monitoring of the companion to detect spectro- 
scopic binaries in several cases), a moderately wide pair in 
a triple system guarantees that the present stellar orbits 
are not disruptive for the planetary system around the iso- 



^" Among 210 members of multiple systems identified by 
Bonavita & Desidera (2007) and Bonavita et al. (|2010p in 
the "Uniform Detectabihty (UD)" sample of Fischer & Valenti 
H2005II . only in 36 cases (18 pairs) were both components of the 
binary system observed. 

^^ 14, if we include HD 97344, which has a close pair of brown 
dwarfs (total mass 0.05 Mq) at a projected separation of about 
2000 AU. 
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late component. A wide stellar triple might also indicate a 
rather unperturbed dynamical history for the system. 

8. Conclusion 

We have presented the discovery of a giant planet with a 
projected mass msmi = 1.49 Mj at 2.6 AU in a moder- 
ately eccentric orbit around HD132563B. This is the first 
planet detected using SARG at TNG as part of the survey 
that monitored about 100 stars that are members of wide 
binaries with similar components. 

The star has a wide companion (HD132563A) at a pro- 
jected separation of 400 AU, which we found to be a spec- 
troscopic binary with a period of about 25-30 yr in a highly 
eccentric orbit. Taking various observational constraints 
into account (unresolved photometry of HD132563A, line 
profile analysis, adaptive optics imaging, astrometric data), 
we concluded that most likely the companion has a mass of 
about 0.5-0.6 Mq. Therefore, HD 132563Bb is one of eight 
planets known to be in stellar triple systems. A tentative 
statistical analysis indicates that the isolated component in 
a hierarchical triple system is not hostile to planet forma- 
tion, with planet frequency as high as for single stars and 
the components of wide binaries. 

HD132563B is a slightly metal-poor star with a mass 
very close to that of Sun. Its companion is less than 200 K 
warmer. Such a small temperature difference, coupled with 
the small convective envelope of the star, makes this sys- 
tem ideal for exploring accretion of metal-rich material of 
planetary origin. The differential abundance analysis (al- 
ready performed in Desidera et al. 12004^ yielded no sig- 
nificant abundance difference, with an error small enough 
to exclude accretion of about 5 M^ of meteoritic material, 
which is the quantity that should have been accreted by 
the Sun during its main sequence lifetime. 

This result, along with the lack of objects in binary 
systems and open clusters showing large enhancements of 
iron abundance that are linked to the evolution of plan- 
etary systems, agrees with the idea that large abundance 
anomalies, comparable to the typical metallicity difference 
between stars with and without giant planets, are rare. 
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Table 6. Abundance differences for binary systems with similar components that host planetary companions. 



System 



Planet 
host 



TeffiA) 

K 



K 



[Fe/H](A) A[Fc/H](A-B) Ref. 



16Cyg 


B 


5745 


62 ±14 


+0.10 


-0.025 + 0.009 


1 


16Cyg 


B 


5765 


39 


+0.05 


0.00 + 0.01 


2 


HD 80606/7 


A 


5700 





+0.46 


-0.01 + 0.11 


3 


HD 80606/7 


A 


5419 


35 


+0.254 


+0.002 + 0.081 


4 


HD 99491/2 


B 


5502 


547 


+0.34 


-0.02 + 0.03 


5" 


HD 99491/2 


B 


5650 


400 


+0.40 


+0.04 + 0.13 


3 


HD 99491/2 


B 


5557 


654 


+0.365 


+0.076 + 0.059 


4 


HD 20781/2 


A 


5774 


518 


-0.06 


+0.05 + 0.03 


6 


HAT-P-1 


B 


6047 


72 


+0.12 


-0.01 + 0.05 


7" 


XO-2 


B 


5500 


160 


+0.47 


+0.02 + 0.03 


ga 


HD 132563 


B 


6168 


184 ± 12 


-0.19 


+0.012 + 0.013 


9 



Refere nces: 1: Laws & Gonzalez (f200T|) : 2: Takeda et al. (f2005|l : 3: Heiter & Luck (f2003)) : 4: Taylor (f2005)) : 5: Valenti & Fischer 
ipnOS)! : 6: Sousa et al. (2008); 7: Bakos et al. (2007) 8; Burke et al. pnD7|) 9: Desidera et al. (2004) 

^■. Valenti & Fischer (2005) actually considered the abundance differences between components of binaries with different temper- 
atures as due to systematic effects in the analysis and introduced a correction to remove them. As these corrections are based 
on the analysis of several pairs, their impact on an individual pair are likely minor. This concern likely applies also to Bakos et 
al. (|2007|l and Burke et al. (|2007[) abundance analysis, which were performed as in Valenti & Fischer (|2005p . 



Table 7. Hierarchical triple systems with planets. The first two columns list the name and the mass of the planet host. 
The 3rd and 4th columns list the individual masses of the two components of the close pair and their semimajor axis or 
projected separation. The 5th column lists the projected separation between the planet host and the distant pair and 
the 6th columns the critical semimajor axis for dynamical stability of planetary companions around the planet host; 
The 7th, 8th and 9th columns list the main planet parameters and finally the 10th column report the references for the 
individual objects. 



Planet Host 



Close Pair 



Wide Pair 



Planet 



Name 



Mass 

Mq 



Comp. Mass 

Mq 



pclose 

AU 



p-wide 

AU 



AU 



Mpi 
Mj 



a.pi 
AU 



Spl 



30 Ari B 
HD 40979 
HD 65216 
HD 132563B 
HD 178911B 
16 Cyg B 
HD 196050 
91 Aqr'^ 



1.16 
1.19 
0.92 
1.01 
1.06 
0.99 
1.15 
1.23 



1.31+0.13*^ 
0.83+0.38 
0.09+0.08 
1.08+0.56^^ 
1.10+0.79 
1.02+0.17 
0.29+0.19 
0.87+0.84 



0.02 
129 

6 
10 

3 

70 
20 
18 



1500 
6400 
253 
400 
640 
850 
511 
2250 



318 
1452 

82 

78 
120 
183 
146 
461 



9.88 
3.83 
1.21 
1.49 
7.35 
1.68 
2.90 
2.90 



1.00 
0.85 
1.47 
2.62 
0.35 
1.68 
2.45 
0.30 



0.29 
0.27 
0.41 
0.22 
0.14 
0.68 
0.23 



Ref 



1,2 

3 

4,5 

6 

7 

8 

5 

9,10,6 



References: 1: Gucnther et al. H2009[l : 2: Morbey & Brosterhus (fT974l; 3: Mugrauer et al. H2007ap . 4: Mugrauer et al. I|2007b|l . 5: 
Eggenberger et al. (2007|; 6: This paper; 7: Tokovinin et al. (2000); 8: Patience et al. flMI} : 9: Raghavan et al. (2006); 10: WDS 
(Mason et al. I200ip 

^■. The critical semimajor axis for dynamical stability was obtained from the projected separation as in Bonavita & Desidera 
l|2007 ). ^■. Minimum mass from spectroscopic orbit. '^: Value of plausible mass, see Sect.O '^■. Planet announced in 2003 (Mitchell et 
al. ED03;) but never published in refereed journals, confirmed at the Conference " Planetary Systems Beyond the Main Sequence" , as 
reported in the Extrasolar Planet Encyclopedia. Stellar mass of 91 Aqr obtained from param interface (Da Silva et al.[2p06) using 
the input parameters from Hekker et al. (|2007ll : stellar masses of the companions from magnitudes in WDS and mass-luminosity 
relation. 
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Table 8. Differential radial velocities of HD 132563 A Table 9. Differential radial velocities of HD 132563 B 



JD -2450000 


RV 


error 


2013.6585 


27.7 


8.3 


2115.4753 


12.3 


8.3 


2327.6477 


18.0 


12.1 


2394.6013 


52.2 


6.8 


2446.5623 


32.5 


10.8 


2774.5525 


71.2 


7.1 


3129.6248 


98.3 


7.1 


3167.5352 


93.8 


6.4 


3574.4756 


119.9 


8.7 


3871.6550 


125.5 


6.1 


3898.5486 


109.5 


5.8 


3901.5701 


116.5 


9.3 


3960.4365 


99.9 


6.3 


3961.4538 


102.5 


6.1 


4099.7730 


113.6 


8.5 


4100.7491 


122.7 


7.0 


4160.6594 


120.7 


7.9 


4161.6741 


101.3 


5.5 


4190.6703 


117.3 


8.0 


4220.6160 


121.4 


6.1 


4221.5886 


129.9 


6.8 


4250.5588 


101.3 


8.8 


4251.5495 


98.6 


8.5 


4252.5885 


87.0 


9.3 


4276.5013 


112.8 


6.8 


4309.4613 


94.8 


8.5 


4311.4002 


68.0 


8.8 


4338.3612 


75.4 


6.6 


4339.3567 


68.1 


6.9 


4369.3521 


112.0 


7.0 


4515.6647 


90.1 


6.7 


4516.6545 


98.9 


11.2 


4544.6446 


74.9 


9.1 


4578.5704 


37.4 


13.1 


4579.5825 


87.9 


8.1 


4605.5027 


66.6 


8.5 


4606.5871 


79.7 


13.0 


4607.5685 


60.9 


9.5 


4609.5616 


69.6 


6.7 


4610.5010 


54.2 


9.1 


4636.5548 


56.5 


8.6 


4669.5064 


65.7 


10.1 


4670.5083 


40.1 


9.2 


4691.4399 


27.5 


6.7 


4849.7753 


-4.9 


9.3 


4961.4343 


-37.1 


5.2 


4962.4406 


-36.8 


7.1 


5046.4187 


-87.6 


10.0 


5054.4395 


-87.8 


8.9 


5449.3704 


-447.4 


12.0 


5600.7471 


-845.1 


7.0 


5643.6802 


-1025.7 


10.4 


5658.6051 


-1061.9 


9.4 



JD -2450000 


RV 


error 


2013.6719 


34.5 


10.8 


2115.4885 


41.4 


11.2 


2327.6601 


21.6 


20.1 


2394.7006 


24.7 


7.9 


2446.5757 


5.1 


12.1 


2774.5654 


-8.4 


7.0 


3129.6370 


35.1 


8.0 


3167.5469 


19.7 


9.6 


3246.4146 


19.9 


10.5 


3371.7370 


26.5 


13.5 


3574.4866 


11.2 


9.9 


3871.6669 


27.7 


6.8 


3898.5605 


-5.6 


5.9 


3901.5821 


-25.9 


9.8 


3960.4494 


-8.1 


6.8 


3961.4655 


-8.3 


7.6 


4099.7852 


-30.3 


12.2 


4100.7608 


-17.8 


9.0 


4160.6721 


-11.6 


9.4 


4161.6960 


-39.2 


8.1 


4190.6825 


-38.1 


9.0 


4220.6272 


-15.3 


8.7 


4221.6003 


-9.8 


6.6 


4250.5700 


-28.6 


10.7 


4251.5607 


3.8 


11.5 


4252.5992 


-21.9 


13.3 


4276.5130 


-30.9 


8.4 


4309.4731 


-19.4 


11.2 


4311.4120 


-33.1 


10.6 


4338.3729 


-18.6 


7.6 


4339.3680 


-20.8 


8.8 


4369.3643 


-7.7 


8.4 


4515.6764 


-6.2 


8.4 


4544.6558 


2.0 


11.3 


4605.5149 


-10.8 


12.0 


4609.5743 


11.8 


8.0 


4669.5191 


29.4 


13.1 


4670.5210 


20.9 


11.1 


4691.4526 


19.9 


9.9 


4961.4455 


30.6 


7.7 


4962.4523 


17.1 


7.4 


5054.4531 


35.8 


13.4 


5434.3939 


18.0 


13.3 


5449.3543 


18.1 


14.5 


5600.7241 


4.9 


7.4 


5643.7032 


-22.2 


9.3 


5658.6217 


-41.3 


11.9 



